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Members  of  the  faculty  who  teach  at  the  undergraduate  and  graduate  levels  and  a number  of 
professional  engineers  and  scientists  whose  primary  activity  is  research  generate  and  conduct  the 
investigations  that  make  up  the  school's  research  program.  The  School  of  Engineering  and  Applied  Science 
of  the  University  of  Virginia  believes  that  research  goes  hand  in  hand  with  teaching.  Early  in  the 
development  of  its  graduate  training  program,  the  School  recognized  that  men  and  women  engaged  in 
research  should  be  as  free  as  possible  of  the  administrative  duties  involved  in  sponsored  research.  In  1959, 
therefore,  the  Research  Laboratories  for  the  Engineering  Sciences  (RLES)  was  established  and  assigned  the 
administrative  responsibility  for  such  research  within  the  School. 

The  director  of  RLES-himself  a faculty  member  and  researcher— maintains  familiarity  with  the 
support  requirements  of  the  research  under  way.  He  is  aided  by  an  Academic  Advisory  Committee  made  up 
of  a faculty  representative  from  each  academic  department  of  the  School.  This  Committee  serves  to  inform 
RLES  of  the  needs  and  perspectives  of  the  research  program. 


In  addition  to  administrative  support,  RLES  is  charged  with  providing  certain  technical  assistance. 
Because  it  is  not  practical  for  each  department  to  become  self-sufficient  in  all  phases  of  the  supporting 
technology  essential  to  present-day  research,  RLES  makes  services  available  through  the  following  support 
groups:  Machine  Shop,  Instrumentation,  Facilities  Services,  Publications  (including  photographic  facilities), 
and  Computer  Terminal  Maintenance. 
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SO.  ABSTRACT  (Cmllnaa  m raaaaaa  al4a  II  naaaaaarr  ma  Uanttfr  hr  black  member) 

This  report  summarizes  the  accomplishments  of  the  final  year  of  a study 
exploring  new  methods  for  balancing,  analyzing  and  designing  flexible  rotating 
shafts.  A technique  for  Identifying  rotor  bearing  parameters  Is  proposed. 

Then  a quadratic  programming  formulation  is  presented  for  Influence  coefficlert 
balancing  with  constraints.  Finally,  a method  Is  developed  for  determining 
the  optimal  axial  location  of  balancing  planes. 
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PROBLEM  STATEMENT  AND  ACCOMPLISHMENTS 


The  goal  of  the  program  is  to  explore  new  techniques  for  balancing, 
analyzing,  and  designing  rotating  shafts.  This  study  began  with  the 
assumptions  that  the  balancing  problem  is  rightfully  classified  as  an 
identification/optimization  problem  and  that  modern  computational  techniques 
hold  great  promise  in  aiding  in  the  solution  of  the  problem.  Through- 
out the  study  emphasis  is  given  to  computational  optimization  formulations. 

The  achievements  of  the  final  year  of  this  study,  which  was  performed 
during  the  period  21  October  1977  - 30  November  1978,  can  be  grouped  Into 
three  categories: 

'Estimation  of  Parameters  in  a Linear  Rotor-Bearing  System  Model 
'The  Determination  of  Rotor  Balance  Weights  by  Quadratic  Programming 
Models 

‘Study  of  the  Sensitivity  of  the  Balanced-Rotor  Residual  Deflections 
to  Deletions  of  Balancing  Weights 
These  efforts  are  summarized  below. 

I.  The  Estimation  of  Parameters  in  a Linear  Rotor-Bearing  System  Model 

For  a system  of  ordinary,  linear  differential  equations  with  constant 
coefficients  which  models  a rotor-bearing  system,  the  problem  of  estimating 
viscous  damping  constants  and  other  parameters  which  appear  in  the  model 
is  approached  from  an  optimization  standpoint.  The  magnitudes  of  the 
constants  are  adjusted  to  give  the  best  agreement  between  some  measured 
values  of  the  system  properties  and  the  values  computed  from  the  model. 

An  estimation  method  is  developed  whereby  a steepest  descent  or  creeping 
random  search  optimization  algorithm  is  applied  to  the  linear  algebraic 
equations  which  govern  the  rotor's  steady  state  response.  The  optimization  is 


conducted.  The  optimization  Is  conducted  at  a single  frequency  and  It  requires 
that  at  least  as  many  displacement  measurements  be  obtained  as  there  are 


unknown  parameters  to  be  estimated.  Because  the  actual  response  characteristics 


of  the  rotor  bearings  are  known  to  be  nonlinear,  the  confinement  of  the 


optimization  to  a single  frequency  should  result  In  a linear  model  which  is 


useful  In  some  range  about  that  frequency.  Of  course.  If  the  linear  model 


is  to  apply  over  a wide  range  of  frequencies,  the  optimization  can  be  per 


formed  using  measurements  taken  at  a number  of  rotor  speeds 


II.  The  Determination  of  Rotor  Balance  Weights  by  Quadratic  Programmini 


In  the  influence  coefficient  method  of  balancing,  the  rotor  response 


is  expressed  as  a linear  function  of  applied  unbalance  weights.  Let  (x,y,z) 


represent  a system  of  coordinates  rotating  at  the  rotor  spin  speed  about 


the  z-axis,  which  is  parallel  to  the  undeformed  shaft  axis.  In  an  isotropic 


system,  the  orbit  of  the  geometric  centroid  of  the  shaft  cross  section  can 


be  considered  to  be  circular  as  it  undergoes  synchronous  whirl  due  to  unbalance 


Let  WX,  WY  represent  the  components  of  deflections  at  various  locations 


along  the  rotor  axis.  Let  UX,  UY  represent  components  of  masses  applied  to 


the  rotor.  Let  [AX],  [AY]  be  matrices  of  influence  coefficients  relating 


the  in-phase  and  out -of -phase  components  of  deflection  to  the  applied 


unbalances  ([AX]  and  [AY]  are  independent  of  the  o 


entatlon  of  the  coordin 


ates) 


Let  WM  b*  i vector  of  displacements  measured  along  the  rotor  which  are 
caused  by  intrinsic  (not  applied)  mass  unbalance.  A typical  rotor 
balancing  procadure  computes  correction  weights  by  minimizing  the  sum 
of  the  squares  of  the  residual  rotor  deflection: 

E • (WM  + W)T(WM  ♦ W)  - (WM  + [A]0)T(WM  + [ A] U) 

In  the  case  in  which  the  elements  of  U are  unconstrained,  the  balance 
weights  are  obtained  by  solving  the  set  of  linear  equations  GRAD(E)  • 0, 

In  which  the  gradient  components  are  derivatives  with  respect  to  the 
elements  of  U.  Often,  the  balance  weights  computed  by  this  method  are  too 
large  to  conveniently  apply  to  the  rotor.  To  overcome  this  problem,  the 
expression  for  E can  be  minimized  subject  to  upper  bounds  on  the  components 
of  UX  and  UY.  Let  X denote  a vector  the  components  of  which  are  sets  of 
pairs  of  real  number  (UX.UY)  for  each  balance  plane.  Carrying  out  the 
multiplications  in  the  expression  for  E results  in: 

E - WO  + C1?  + XT(D]X 

WO  is  the  sum  of  the  squares  of  the  measured  deflections. 

C is  the  cost  vector  (a  function  of  WM  and  [A]). 

[D]  is  the  quadratic  matrix  which  is  assumed  to  be  positive  definite  ([D]  is 
a function  of  [A]). 

Let  UM  be  a vector  of  the  upper  and  lover  bounds  on  the  elements  of  X. 

- UM  < X < UM 

By  lntroductlng  a shift  of  coordinates,  these  constraints  can  be  converted 
to  the  fora  0 \ ' <^UM  ♦ UM  In  which  • X + UM. 

In  this  study,  a quadratic  programming  routine  was  employed  to 
demonstrate  the  usefulness  of  the  quadratic  programing  algorithm  to  rotor 
balancing. 
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ly  attractive  exact  method  of  finding  the  correct  answer  does  not  appear 
to  be  possible.  Suppose  that  the  optimal  choice  of  n of  N possible  planes 


N' 

(n<N)  is  sought;  then,  there  will  be  nV'(^_n)  ; possible  choices.  As  N 
increases,  an  exhaustive  search  rapidly  becomes  impractical  (e.g.  for 
N*20,  there  are  184,756  possible  choices  of  10  planes).  A random  search 
algorithm  could  be  a feasible  approach;  for  example,  the  probability  of 
finding  a solution  in  the  best  1%  of  all  solutions  is  better  than  99% 
with  a random  sample  of  500  trials.  However,  running  even  500  problems  is 
a fairly  unattractive  alternative. 

Because  only  n-plane  solutions  which  do  not  too  greatly  increase 
the  rotor  residual  deflection  from  the  N-plane  level  are  of  interest,  a 
region  near  the  N-plane  optimum  may  exist  where  variable  interactions 
can  be  neglected  in  estimating  the  effects  of  balance  plane  deletions.  Any 
technique  which  attempts  to  consider  the  interactions  of  delected  balance 
planes  would  entail  examination  of  all  or  a great  number  of  the  N!/(N-n)!n! 
possibilities.  By  considering  only  independent  effects,  a small 
fraction  of  the  great  number  of  possibilities  can  be  segregated  for  exact 
solution. 

To  estimate  the  effects  of  balance  plane  delectlons,  two  studies  were 
conducted.  In  the  first,  the  unconstrained  solution  with  N-planes  was 
obtained;  then  N additional  problems  were  solved  in  which  the  upper  bounds 


on  on*  of  the  N-planes  was  reduced  either  to  zero  or  to  a very  small 
value.  From  the  change  in  the  objective  function  and  the  value  of  the 
Lagrange  multiplier  for  each  deleted  variable,  the  most  insensitive 
balance  plane  can  be  estimated.  These  estimated  selections  are  then 
compared  with  actual  n-plane  solutions  determined  by  exhaustive  or  random 
searches. 

In  the  second  study,  the  constraint  limits  were  allowed  to  approach 
zero  simultaneously,  and  the  order  in  which  constraints  become  active, 
the  change  in  the  objective  function  as  constraints  decreased,  and  the 
changes  in  the  Lagrange  multipliers  were  studied.  This  Information  ms 
used  to  estimate  balance  plane  sensitivity, as  in  the  first  study.  These 
results  indicate  that  independent  effects  can  be  used  to  estimate  the 
sensitivity  of  the  rotor  balance  level  to  plane  delections. 
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The  University  of  Virginia's  School  of  Engineering  and  Applied  Science  has  an  undergraduate 
enrollment  of  approximately  1,300  students  with  a graduate  enrollment  of  approximately  500,  There  are 
125  faculty  members,  a majority  of  whom  conduct  research  in  addition  to  teaching. 

Research  is  an  integral  part  of  the  educational  program  and  interests  parallel  academic  specialties. 
These  range  from  the  classical  engineering  departments  of  Chemical,  Civil,  Electrical,  and  Mechanical  and 
Aerospace  to  departments  of  Biomedical  Engineering,  Engineering  Science  and  Systems.  Materials 
Science,  Nuclear  Engineering  and  Engineering  Physics,  and  Applied  Mathematics  and  Computer  Science. 
In  addition  to  these  departments,  there  are  interdepartmental  groups  in  the  areas  of  Automatic  Controls  and 
Applied  Mechanics.  All  departments  offer  the  doctorate;  the  Biomedical  and  Materials  Science 
Departments  grant  only  graduate  degrees 

The  School  of  Engineering  and  Applied  Science  is  an  integral  part  of  the  University  (approximately 
1,530  full-time  faculty  with  a total  enrollment  of  about  16,000  full-time  students),  which  also  has 
professional  schools  of  Architecture,  Law,  Medicine,  Commerce,  and  Business  Administration.  In  addition, 
the  College  of  Arts  and  Sciences  houses  departments  of  Mathematics.  Physics,  Chemistry  and  others 
relevant  to  the  engineering  research  program.  This  University  community  provides  opportunities  for 
interdisciplinary  work  in  pursuit  of  the  basic  goals  of  education,  research,  and  public  service. 


